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Abstract 

A  ternary  phase-field  model  was  developed  that  is  linked  directly  to  commercial 
CALPHAD  software  to  provide  quantitative  thermodynamic  driving  forces.  A  re¬ 
cently  available  diffusion  mobility  database  for  ordered  phases  is  also  implemented 
to  give  a  better  description  of  the  diffusion  behavior  in  alloys.  Because  the  targeted 
application  of  this  model  is  the  study  of  precipitation  in  Ni-based  superalloys,  a 
Ni-Al-Cr  model  alloy  was  constructed.  A  detailed  description  of  this  model  is  given 
in  the  paper.  We  have  considered  the  misfit  effects  of  the  partitioning  of  the  two 
solute  elements.  Transformation  rules  of  the  dual  representation  of  the  7  +  7'  mi¬ 
crostructure  by  CALPHAD  and  by  the  Phase-Field  are  established  and  the  link  with 
commercial  CALPHAD  software  is  described.  Proof-of-concept  tests  were  performed 
to  evaluate  the  model  and  the  results  demonstrate  that  the  model  can  accurately 
describe  the  compositional  and  structural  changes  associated  with  precipitation  of 
the  7'  phase.  Uphill  diffusion  of  A1  is  observed  in  a  few  diffusion  couples  showing 
the  significant  influence  of  Cr  on  the  chemical  potential  of  Al.  Possible  applications 
of  this  model  are  discussed. 

Key  words:  Ni-Al-Cr  alloys,  ternary  phase-field  models,  precipitation,  commercial 
CALPHAD  software,  uphill  diffusion 
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1  Introduction 


Traditionally,  precipitate  processes  have  been  studied  in  a  largely  empirical 
fashion,  with  samples  being  prepared,  carefully  heat  treated,  and  character¬ 
ized.  In  more  recent  years,  there  has  been  an  explosion  of  materials  modeling 
techniques  whose  predictive  capabilities  can  be  a  guide  for  experimental  efforts 
to  increase  the  yield  above  that  for  simple  trial-and-error  methods.  Therefore, 
it  would  be  desirable  to  have  an  engineering  tool  that  could  predict  quantita¬ 
tively  the  microstructural  evolution  of  an  alloy  under  specified  thermodynamic 
conditions  during  heat  treatment. 

The  phase-field  method  has  shown  great  potential  for  predicting  realistic  mi¬ 
crostructures.  However,  most  phase-field  simulations  performed  to  date  have 
been  limited  to  binary  systems  such  as  Ni-Al  [1,  2,  3,  4,  5,  6]  and  most  of 
the  models  employed  were  qualitative.  Recently,  attempts  have  been  made  to 
develop  quantitative  phase-field  models  by  specifying  the  fundamental  param¬ 
eters  using  thermodynamic  and  diffusivity  databases  [7,  8,  9,  3,  4,  10,  6].  Since 
Ni-based  superalloys  contain  more  than  a  dozen  elements,  the  development  of 
quantitative  multi-component  phase-field  models  is  crucial. 

The  Ni-Al-Cr  system  is  the  most  important  ternary  alloy  for  Ni-based  su¬ 
peralloys,  and  it  provides  the  basis  for  simplifying  the  complex  chemistry  of 
more  complicated  Ni-based  superalloys  [11].  The  importance  of  this  system  is 
reflected  by  the  numerous  experimental  studies  on  its  phase  equilibrium  and 
thermodynamics;  for  a  general  review,  see  [11,  12],  Based  upon  these  extensive 
experimental  studies,  critically  assessed  thermodynamic  and  kinetic  mobility 
databases  are  now  available  in  the  literature  [13,  14,  11,  15]. 

In  this  paper  we  present  a  ternary  phase-field  model  that  links  directly  to 
commercial  CALPHAD  software  to  obtain  the  thermodynamic  driving  forces, 
and  to  a  recently  developed  diffusion  mobility  database  [16]  that  describes  the 
contribution  from  the  disordered  fee  phase  as  well  as  the  ordered  Z/l2  phase. 
Description  of  the  model  is  presented  in  Sec.  2  followed  by  the  results  of  some 
proof-of-concept  simualtions;  a  brief  discussion  of  the  capability  of  the  model 
in  Sec.  3. 


2  The  ternary  phase-field  method 


The  ternary  phase-field  model  is  linked  directly  to  commercial  CALPHAD 
software.  We  have  implemented  a  diffusion  mobility  database  that  became 
available  very  recently.  This  database  allows  us  to  have  a  better  description 
of  diffusion  behavior  near  the  7/7'  interfaces  and  inside  the  ordered  7'  phase. 
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In  the  phase-field  formulation,  we  have  also  considered  the  misfit  effect  due  to 
the  partitioning  of  the  two  solute  elements. 


2.1  Dual  representation  of  the  7  +  7'  microstructure  by  CALPHAD  and 
Phase-Field 


There  is  an  incompatibility  between  the  descriptions  of  ordered  domain  struc¬ 
ture  in  the  CALPHAD  and  Phase-Field  methods.  After  reviewing  each  method, 
we  will  develop  the  transformation  rules  from  one  system  to  the  other  that 
is  necessary  for  our  link  between  our  Phase-Field  model  and  the  commercial 
software  based  on  CALPHAD  method  that  we  employ. 


2.1.1  CALPHAD  representation  of  the  7  +  7'  microstructure 

In  the  CALPHAD  method,  sublattices  are  used  to  describe  ordering  that  can 
not  be  described  by  the  regular  solution  [17].  In  the  regular  solution  approx¬ 
imation,  all  species  involved  are  assumed  to  be  sufficiently  similar  in  terms 
of  size,  shape,  electronegativety,  etc.,  that  random  mixing  of  the  solutes  is  a 
reasonable  approximation.  If  the  species  differ  significantly,  random  mixing  is 
not  a  good  approximation  and  as  a  consequence  the  solution  can  not  be  well 
described  by  the  ideal  entropy  of  mixing.  In  this  case,  different  species  may 
occupy  different  sublattices  (a  regular  solution  is  usually  assumed  within  each 
sublattice)  and  the  distribution  of  the  species  are  described  by  their  site  frac¬ 
tions  in  the  sublattices.  Depending  on  the  specific  components  in  a  system,  the 
number  of  sublattices  needed  for  an  appropriate  description  of  the  system  are 
usually  different.  For  the  7  +  7'  microstructure  in  Ni-based  alloys,  descriptions 
using  two  or  four  sublattices  can  be  found  in  the  literature  [18,  19,  14]. 

An  example  of  a  binary  model  using  four  sublattices  to  describe  the  Ni-Al 
system  is  given  below.  A1  atoms  occupy  the  corners  of  the  fee  lattice  in  one  of 
the  sublattices  while  the  Ni  atoms  occupy  the  center  (left  side  of  Fig.  1(a)). 
There  are  three  variants  (or  domains)  of  this  sublattice;  one  is  illustrated  on 
the  right  side  of  Fig.  1(a)  where  the  two  A1  atoms  are  on  the  opposing  face 
centers  along  the  [100]  direction  and  the  Ni  atoms  occupy  the  rest  of  the 
lattice  .  The  two  A1  atoms  in  other  two  variants  similarly  occupy  opposing 
face  centers  but  along  the  [010]  and  [001]  directions,  respectively.  Projections 
of  the  distribution  of  atoms  along  [001]  direction  are  illustrated  in  Fig.  1(b) 
for  the  four  different  ordered  domains. 

The  site  fractions  for  the  sublattices  are  subject  to  the  normalization  condition 
that  summing  the  fractional  occupations  over  all  the  sublattices  must  yield 
the  concentration  field: 
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(1) 


7 Y,yS  M)  =  x(r,t) 

^s= 1 

Following  the  discussion  of  the  binary  system,  there  are  four  different  domains 
in  terms  of  each  solute  elements  depending  on  how  they  are  partitioned  in 
the  fee  structure.  The  site  fractions  for  the  ternary  alloys  are  subject  to  the 
normalization  condition: 

7E^(bf)  =  s«(bi)  (2) 

S—  1 

It  should  be  pointed  out  that  the  partitioning  of  the  two  solute  elements  must 
be  matched  with  each  other  in  the  precipitate  to  maintain  the  Ll2  structure. 
If  the  sublattices  of  Cr  and  A1  were  different,  the  Al/Cr  ratio  would  not  be 
the  same  for  all  sublattices  and  an  additional  ordering  reaction  would  occur 
resulting  from  this  loss  of  symmetry.  For  example,  Fig.  2  shows  what  would 
happen  if  one  of  the  base  centered  lattice  site  was  not  equivalent  to  the  two 
face  centered  lattice  sites  with  respect  to  Cr  concentration.  The  difference  in 
Cr  concentration  in  the  base  centered  site  would  make  this  site  distinct  from 
the  face  centered  sites  and  a  tetragonal  crystal  structure  Ll2  would  result. 


2.1.2  Phase-fi,eld  representation  of  the  7  +  7'  micro  structure 

In  order  to  describe  the  binary  Ni-Al  phase-field  model,  we  require  three  long- 
range  order  ( Iro )  parameters  rf(r,t)  with  p  =  1,2,3  (where  r  and  t  are  the 
spatial  coordinates  and  time)  and  a  concentration  field  x(r ,  t).  These  variables 
fully  describe  the  partitioning  of  A1  on  the  four  sublattices  of  the  Ll2  struc¬ 
ture  (see  Fig.  1)  and  hence  the  composition  inhomogeneity  in  the  two-phase 
microstructure  (7  +  V)  t^- 

Let  p  (r,  t )  denote  the  probability  of  finding  solute  atom  at  site  r  at  time  t ,  the 
order  parameters  rf{r,t)  are  defined  as  Fourier  coefficients  in  an  expansion  of 
the  probability  p(r,t)  as  in  Eq.  3: 


P  (r,  t) 

k1 


=x  (r,  t)  jl  +  rj1  (r,  t )  exp  [jik1  •  rj 

+?f  (r,  t )  exp  [ik2  •  rj  +  ?/3  (r,  t )  exp  [ik3  •  rj  j 

=  —  [100] .  k2  =  —  [010] ,  k3  =  —  [001] 
a  a  a 


(3) 


The  concentration  field  x{r ,t)  can  then  be  identified  as  the  probability  of 


4 


finding  solute  atom  at  site  r  at  time  t  in  the  completely  random  alloy,  as  in 
Eq.  4: 


V1  (r,  t)  =  V2  (r,  t)  =  V3  (r,  t)  =  0  =>  p  (r,  t)  =  x  (r,  t)  (4) 

Let  rj°  denote  the  equilibrium  Iro  parameter  at  a  given  composition  and  tem¬ 
perature.  Then  four  energetically  equivalent  antiphase  domains  of  the  ordered 
Ll2  structure  [20]  can  be  defined  by  specifying  the  order  parameters  (771,  rj2,  r\ 3) 
as  the  constant  values  in  Eq.  5: 


domain  rf  (r,  t)  rf  (r,  t)  rf  (r,  t) 


1  T]° 

2  -rj° 

3  ~7]° 

4  n° 


f]°  T]° 

—rj°  rf 

rf  —rj° 

—rj°  — rf 


(5) 


For  example,  the  first  domain  is  characterized  by  enriched  solute  at  the  eight 
corners  of  the  cubic  cell  and  depleted  solute  at  the  six  face  centers  (Fig.  1(a)), 
as  in  Eq.  6.  The  other  three  domains  have  enriched  solute  at  the  two  opposing 
fee  face  centers  normal  to  the  [001],  [010],  and  [100]  directions  respectively 
and  depleted  solute  at  the  other  face  centers  and  corners  of  the  cubic  cell. 


dM) 


t?2(M)  =  ??3(r,t)  =  V°  =>  (6) 

x(t,  t )  [1  +  3^°] ,  k1  •  r  =  k2  •  r  =  k3  •  r  =  2mr 

x(t ,  t)  [1  —  770] ,  k1  •  r  =  k2  •  r  =  (2 n  +  l)7r,  k3  •  r  =  2nn 

< 

x(r.  t )  [1  —  770] ,  k2  •  r  =  k3  •  r  =  (2 n  +  l)7r,  k1  •  r  =  2rm 

x(r.  t)  [1  —  rj°] ,  k3  •  r  =  k1  •  r  =  (2 n  +  l)7r,  k2  •  r  =  2n7r 


To  extend  this  binary  analysis  to  the  ternary  Ni-Al-Cr  alloys,  an  additional  set 
of  four  field  variables  are  needed  to  describe  the  partitioning  of  Cr  in  the  Ll2 
structure  and  the  composition  inhomogeneity  of  Cr  in  the  two  phases.  Conse¬ 
quently,  we  need  a  total  of  eight  field  variables  to  describe  the  microstructure 
in  a  ternary  phase- field  model,  i.e.  six  Iro  parameters  r]f(r,t)  with  p  =  1,2,3 
and  i  =  1,2;  here  i  =  1  represents  A1  and  i  =  2  represents  Cr.  Similarly  we 
need  two  concentration  fields  x,(r,t)  to  describe  the  two  solutes  A1  and  Cr. 
The  probabilities  of  finding  an  A1  or  Cr  solute  atom  are  now  extended  by 
analogy  for  ternary  case  as  in  Eq.  7 
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(7) 


Pn  (r,  t)  =  xn  (r,  t)  {l  +  77*  (r,  t)  exp  [ik1  •  r] 

Wn  (r,  t)  exp  [?:k2  •  r]  +  rfn  (r,  t)  exp  [ik3  •  r]  } 


2.1.3  Transformation  of  the  Dual  Representation  between  CALPHAD  and 
Phase-Field 

In  the  phase-field  representation  of  the  7+7'  microstructure  in  binary  systems, 
we  used  the  variables  771  (r, t) ,  rj2{r,t)1  if{r,t)1  and  x(r ,t).  In  the  CALPHAD 
represenation,  the  free  energy  is  parameterized  by  the  site  fractions  of  the 
solute,  y 1  (r,  t).  y 2  (r,  t).  y 3  (r,  f),  and  y4  (r,  t).  These  are  completely  equivalent 
formulations  that  describe  ordered  alloys.  While  the  order  parameters  have  an 
explicit  meaning  for  the  degree  of  ordering  in  the  alloy,  the  site  fractions  have 
a  clear  correspondence  for  the  occupation  of  the  various  sublattices. 

We  must  transform  between  these  equivalent  formulations  in  order  to  employ 
the  free  energies  computed  with  the  thermodynamic  databases  in  the  phase- 
field  simulations.  The  correspondence  between  site  fractions  of  the  solute  ys 
and  the  phase-field  parameters  x  and  rf  is  shown  in  Eq.  8  where  each  domain 
ys  represents  the  site  fractions  of  the  solute  A1  on  sublattice  s  [3]. 


y 1 M) 
y2  (r,  t) 
y 3  (r,  t) 
y 4  M) 


x  (r,  t) 

1  +  7 g  (r,  t)  +  rf  (r,  t)  +  rf  (r,  t) 

x  (r,  t) 

1  -  V1  (r,  t)  -  if  (r,  t)  +  rf  (r,  t) 

x  (r,  t) 

1  -  rf  (r,  t)  +  rf  (r,  t)  -  rf  (r,  t) 

x  (r,  t) 

1  +  V1  (r ,t)  -  rf  (r ,t)  -  rf  (r,i)j 

(8) 


Transformation  between  the  dual  representations  for  ternary  systems  can  be 
obtained  by  a  straightforward  extention  of  Eq.  8: 


Vn  (r,  t)  =  xn  (r,  t) 
IJn  Ob  t)  =  xn  (r,  t) 

lli  (r,  t)  =  xn  (r,  t) 

Vn  (r’  *)  =  xn  (r,  t) 


1  +  Vn  (r’ t)  +  rji  (r,  t)  +  rfn  (r,  t)] 

1  -  Vn  (r’  f)  ~  Vn  (r’ t)  +  Vn  (r^  0 

1  -  Vn  (r^  t)  +  r)l  (r,  t)  -  r/3  (r,  t) 

1  +  Vn  (r^  t)  -  rg  (r,  t)  -  rg  (r,  t) 


(9) 
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2.2  Ternary  Phase-Field  Evolution  Equations 


The  microstructural  evolution  can  be  simulated  by  solving  the  time-dependent 
Ginzburg-Landau  equation  for  the  six  Iro  parameters  using  the  non-linear 
Cahn-Hilliard  diffusion  equation  for  the  two  concentration  fields: 


drji 

dt 

dxj 

dt. 


y2V 

m 


Y.M» 


(10) 

(11) 


Here  L  is  the  structural  relaxation  constant,  are  the  diffusion  mobilities, 
Vm  is  the  molar  volume,  and  F  is  the  total  free  energy  of  the  system. 

For  the  coherent  precipitation  under  consideration,  the  total  free  energy  con¬ 
sists  of  the  elastic  strain  energy  (Fei)  and  the  chemical  free  energy  (Fch): 


F  —  Fet  +  Fch 


(12) 


According  to  gradient  thermodynamics  [21],  the  chemical  free  energy  should 
consist  of  the  local  chemical  free  energy  and  the  gradient  terms  arising  from 
compositional  and  structural  inhomogeneities,  i.e., 


Fch=  f  f(xAi,xCr,rfiu,rfcr,T)  +  ]-  ^  UXij  VxiVxj)  (13) 

V  ^  a  a—  A  i 


i,j-Al,Cr 


+\Kn  E  E  (yrfiy 

Z  i=Al,Cr  P 


dV 


where  /  (x^i,  xcr ,  Vai-i  Vcn  T)  is  the  local  chemical  free  energy  that  defines  the 
basic  thermodynamic  properties  of  the  system,  nXij  and  are  gradient  energy 
coefficients,  and  it  is  understood  that  p  =  1,2,3  for  the  three  components  of 
the  Iro  parameters. 
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2.3  Elastic  Strain  Energy 


The  elastic  strain  energy  arises  because  of  the  lattice  mismatch  between  the 
7  and  7'  phases.  The  misfit  strain  is  assumed  to  depend  only  on  the  compo¬ 
sitions,  and  the  linear  terms  in  the  expansion  of  the  strain  are: 


eii(r)  =  12AXP  W  6ij(P) 

P 

(14) 

A xp  (r)  =  xp  (r)  -  x°p 

(16) 

Here  p  =  1  implies  A1  and  p  =  2  implies  Cr.  The  quantity  Axp  measures  the 
deviation  of  the  local  concentration  of  a  solute  from  its  total  site  fraction  x°, 
and  e°j  (p)  is  the  Stress  Free  Transformation  Strain  (SFTS)  tensor.  For  the 
7  -*  7'  transformation,  the  SFTS  is  characterized  by  a  volume  change  that 
can  be  approximated  by  Vegard’s  law: 


e°j  (r)  =  { 


tp...i=  3 

0  ...i  +  j 

da  (xaii  •Ect) 
a0dx  „ 


(16) 


Here  ep  is  the  concentration  coefficients  of  the  crystal  lattice  parameter,  a(xAi ,  %cr) 
is  the  lattice  parameter  of  a  solid  solution  with  concentration  xai  and  xcr,  and 
a0  is  the  lattice  parameter  of  pure  solvent.  The  elastic  energy  of  a  two-phase 
mixture  can  then  be  calculated  as  a  function  of  the  concentration  profiles: 


V  _  *  _ 

Eel  V ^  1  ^kl  ip)  A^  ( P )  (17) 

1  p=  1 

+  \cijkl  E  E  eh  (p)  eh  (?) ip)  Ax  (q) 

Z  p=lg=l 

-JEE  /  t^-T3  bpq  (n)  iAx  (p)Ys{Ax  («)}g 
Z  p=lq=lJ  (2?r) 

Here  (A.)  is  the  average  volume  of  (A).  V  is  the  total  volume  of  the  system, 
Cijki  is  the  elastic  moduli  tensor,  the  notation  f'  indicates  that  the  point 
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n  =  0  is  excluded  from  the  integration,  and  Bpq  (n)  is  a  two-body  interaction 
potential  given  by: 


Bpq{ n)  =  riia-j  ( p )  Qjk  (n)  aki{q)ni  (18) 

Here  n  =  is  a  unit  vector  in  reciprocal  space  and  n,  is  the  ith  component, 
a°j  =  Cijkitkl  ( p )  and  Qt]  (n)  is  a  Green  function  tensor  which  is  inverse  to  the 
tensor  0  A  (n)  =  Cijkinkni,  {Ax  (p)}s  is  the  Fourier  transform  of  Ax  ( p ),  and 
{Ax  (p)}*  is  the  complex  conjugate  of  {Ax  (p)}g. 

2.4  Linking  the  chemical  free  energy  from  thermodynamic  databases  to  phase- 
fi,eld  simulations 


In  the  CALPHAD  method,  the  local  chemical  free  energy  is  expressed  as  a 
sum  of  a  reference  free  energy,  an  ideal  mixing  free  energy,  and  excess  mixing 
free  energy: 


G  (r,  T)  =  Gref  (r,  T)  +  Gid  (r,  T)  +  Gex  (r,  T)  (19) 

The  reference  free  energy  for  substitutionally  disordered  solid  solutions  is: 


Gref  (r,  T)  =  J2  Vi  (r)  Vj  (r)  vl  (r)  Vi  (r)  G°ijki  (T)  (20) 

i— 1 j= 1 k= 1 1=1 

Here  N  is  the  total  number  of  elements  in  the  alloy,  and  yf  (r)  is  the  site 
fraction  of  element  n  on  sublattice  s  with  normalizations  yp  (r)  =  1  and 
4'Z)s=i  yp  (r)  =xp ;  that  is,  summing  the  fractional  occupancies  over  all  the  al¬ 
loying  elements  on  each  sublattice  must  equal  one,  and  summing  the  fractional 
occupancies  over  all  sublattices  for  each  alloying  element  must  equal  the  mole 
fraction  of  that  element.  The  functions  G°]kl  ( T )  are  empirically  fitted  pa¬ 
rameters  tabulated  as  functions  of  temperature  in  thermodynamic  databases. 
The  parameters  themselves  represent  the  free  energies  associated  with  those 
phases  in  which  each  sublattice  is  occupied  by  a  single  element,  the  so-called 
’end-member  compounds’.  These  free  energies  are  mixed  according  to  the  site 
fractions  to  describe  a  particular  alloy.  The  ideal  mixing  free  energy  Gld  (r,  T) 
is: 


Gu  (r,  T)  =  RT  £  /■  ■£  y‘n  (r)  In  K  M] 


s=l  n= 1 


(21) 
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Here  R  is  the  gas  constant,  fs  is  the  site  fraction  of  sublattice  s.  This  ex¬ 
pression  gives  the  entropy  of  mixing  of  an  ideal  solution,  enumerating  all  the 
equally  possible  rearrangements  of  non-interacting  atoms  within  each  sublat¬ 
tice.  Finally,  the  excess  free  energy  of  mixing  Gex  (r,  T)  is: 


G“  (r,  T)  =  £  £  £  y;  (r)  •£  y|  (r)  y*  (r)  (22) 

s—lt^sp—l  q— 1 

Nipq 

■  V  K  (r)  -  y-„  (r)  PG-  (m.  T) 

m=  1 


Here  Gs^q  (m,  T)  are  empirically  fitted  parameters  tabulated  as  functions  of 
temperature  in  thermodynamic  databases. 


The  local  chemical  free  energy  in  Eq.  19  is  expressed  as  a  function  of  the 
site  fractions  of  the  elements  on  each  of  the  four  sublattices.  However,  in  the 
phase-field  description  the  local  chemical  free  energy  is  a  function  of  local 
concentrations  and  Iro  parameters.  The  two  set  of  parameters  are  related  to 
each  other  through  Eq.  9. 


We  have  chosen  to  link  our  phase-field  code  to  the  commercial  CALPHAD 
software  Pandat  (CompuTherm,  LLC).  The  coupling  between  the  phase- field 
simulation  and  the  commercial  thermodynamic  database  involves  the  evalua¬ 
tion  of  the  partial  derivatives  of  the  local  chemical  free  energy  with  respect  to 
the  eight  site  fractions  (ysn  (r))  by  Pandat.  These  partial  derivatives  are  then 
used  to  evaluate  the  partial  derivatives  of  the  local  chemical  free  energy  with 
respect  to  the  phase-field  parameters,  i.e.  concentration  parameters  (xn  (r)) 
and  the  Iro  parameters  ( 77 d  (r) ,  rfa  (r) ,  rfa  (r)),  through  the  chain  rules: 


df_ 

dxp 


4  2 


EE 


dG  dysn 
d'!Jn  dxp  ' 


..p=  1,2,3 


df 


4  2 


EE 


dG  dysn 
dysn  drj%  ' ' 


.  k  =  1,2,3 


(23) 


2.5  Diffusion  Mobility 


The  diffusion  mobility  varies  with  composition  and  structure.  Due  to  very 
limited  data  on  the  diffusion  mobility  from  the  ordered  phase,  most  simula¬ 
tions  to  date  have  used  the  diffusion  mobility  of  the  disordered  fee  phase  to 
represent  that  of  the  overall  behavior  of  the  alloys.  In  other  words,  the  contri¬ 
bution  from  the  ordered  phase  is  totally  ignored  although  it  is  well  understood 
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that  the  diffusion  behavior  and  therefore  the  closely  related  diffusion  mobility 
in  the  ordered  phase  are  quite  different  from  that  in  the  disordered  fee  phase. 
Most  recently,  Campbell  has  developed  a  diffusion  mobility  database  that  is 
capable  of  describing  the  diffusion  mobilities  in  the  7'  phase  in  the  Ni-Al-Cr 
system  [16].  In  this  work,  we  adopted  this  database  for  the  description  of  the 
diffusion  mobilities  in  the  Ni-Al-Cr  system. 

The  diffusion  mobility  is  expressed  as  [22] 

MaIAI  —  TZ~xAl  \xCrxAlBcr  +  (1  —  xAl )2  Bai  +  xAlxNiB]sri\  (24) 

*m 

MaiCt  —  Tl~xAlxCr  [—  (1  —  xAl)  Bai  ~  (1  —  xCr)  Bcr  + 

* m 

McrCr  —  7~p  xCr  \xAIxCtBaI  +  (1  _  xCr)2  Bcr  +  xNixCrBNi ] 

where  R  is  the  atomic  mobility  of  species  i.  The  atomic  mobility  is  usually 
expressed  as: 

B,  =  ^exp(-AQ;/RT)  (25) 

Here  R  is  the  gas  constant  and  Q*  is  the  total  activation  energy  of  species  i 
that  consists  of  contributions  from  the  disordered  (A  Qfs)  and  ordered  (A  Q°rd) 
states  as  described  below, 

A  Q*  =  A  Qfs  +  AQ°rd  (26) 

The  contribution  from  the  disordered  state,  A Qf s,  and  its  composition  depen¬ 
dence  can  be  obtained  based  on  Engstrom  and  Agren’s  work  [15].  The  contri¬ 
bution  from  the  ordered  state  is  described  below  as  in  Campbell’s  work  [16], 

*QZd  =  E  E  ^QZMvi  ~  x<*i}  +  £  £  E  ^QTrMv^l  - 

i  i  j  k 

+  E  AQ™ij  [yf  y?  y%  -  xixJxk ] ,  (27) 

i  j  k 

where  the  index  m,  i,j,k  =  1,  2,  3  representing  the  three  components  Al ,  Cr,  Ni 
in  the  ternay  system,  yf  is  the  site  fraction  of  component  i  on  the  a  sublattice 
that  can  be  readily  converted  from  the  current  4-sublattice  formulation  [19]. 
A A Q™.k,  are  the  contributions  to  the  activation  energy  for  compo¬ 

nent  m  as  a  result  of  the  chemical  ordering  of  the  i.  j,  and  k  atoms  on  the  two 
sublattices.  These  parameters  are  well  documented  in  Campbell’s  work  [16]. 


3  Results  and  Discussion 


We  have  presented  a  ternary  phase-field  model  that  links  directly  to  commeri- 
cal  CALPHAD  software  and  uses  a  comprehensive  diffusion  mobility  database 
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that  describes  the  contribution  from  the  disordered  fee  phase  as  well  as  the 
ordered  Ll2  phase.  Proof  of  principal  simulations  were  performed  to  determine 
if  the  model  can  accurately  describe  the  compositional  and  structural  changes 
associated  with  the  precipitation  behavior  in  7  — >-  7'  phase  transformation. 
The  uphill  diffusion  behavior  in  Ni-Al-Cr  ternary  systems  is  simualted,  which 
is  a  typical  phenomenon  of  systems  with  more  than  two  components.  In  the 
end  of  this  section,  some  further  opportunities  for  applying  this  model  are 
discussed. 

As  mentioned  in  Sec.  2.1.2,  the  domain  structure  of  the  two  solute  elements 
must  match  to  maintain  the  desired  Ll2  structure  in  the  precipitate.  This 
means  that  the  profiles  of  three  Iro  parameters  for  Al  need  to  overlap  with 
those  for  Cr.  To  examine  if  the  model  can  properly  describe  precipitation, 
simulations  were  run  to  observe  the  ordering  during  the  growth  of  a  single 
precipitate  from  a  small  seed  in  a  supersaturated  matrix.  In  these  simulations, 
we  assume  equilibrium  compositions  for  the  seed;  also  it  is  fully  ordered  with 
Vn  =  Vn  ~  Vn  =  V°  where  n  =  1,2  represent  respectively  solute  Al  and 
Cr.  Two  different  bulk  compositions  in  the  matrix  were  chosen  to  investigate 
the  possible  effect  of  varying  the  concentration  partitioning  on  the  ordering. 
The  compositions  were  Ni-lOAl-lOCr  at.%  and  Ni-15Al-5Cr  at.%,  and  the 
temperature  was  873K. 

Fig.  3  shows  the  growth  of  the  precipitate  and  the  evolution  of  the  concentra¬ 
tions  and  the  Iro  parameters  77*  forAi  ( n  =  1)  and  Cr  ( n  =  2)  in  Ni-lOAl-lOCr 
at.%.  It  was  observed  that  7/ ?  and  follow  the  evolution  of  77*  very  closely. 
The  three  representative  snapshots  clearly  show  that  the  two  order  parame¬ 
ters  r}\  and  772  follow  one  another.  The  compositions  in  the  precipitate  show 
an  enrichment  in  Al  and  a  depletion  in  Cr  as  expected  from  the  ternary  phase 
diagram.  This  demonstrates  that  the  model  can  adequately  describe  the  com¬ 
positional  and  structural  changes  associated  with  precipitation. 

For  alloy  Ni-15Al-5Cr  at.%,  the  partitioning  of  Cr  in  the  two  phases  7  and 
7'  is  about  equal  based  on  the  ternary  phase  diagram.  Therefore,  the  growth 
of  the  7'  precipitate  is  expected  to  involve  mainly  diffusion  of  Al  while  the 
diffusion  of  Cr  should  be  insignificant.  This  is  indeed  the  case  in  the  results  of 
the  simulation  as  shown  in  Fig.  4.  The  Cr  composition  is  largely  the  same  in 
the  two  phases  except  near  the  interfaces  where  significant  variation  is  present. 
However,  the  lack  of  Cr  partitioning  does  not  affect  the  ordering  behavior  at  all 
and  the  two  order  parameters  still  closely  follow  each  other  with  the  growth 
of  the  precipitate.  This  is  a  further  indication  that  the  model  can  properly 
describe  the  Ll2  structural  change  associated  with  the  precipitation  process. 

Simulation  tests  were  extended  to  include  multiple  precipitate  particles,  where 
nuclei  were  introduced  into  individual  cells  randomly,  with  the  probability  re¬ 
lated  to  the  corresponding  nucleation  rate  of  those  cells  [23].  The  order  pa- 
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rameters  for  each  of  the  nuclei  are  also  assigned  randomly  to  maintain  roughly 
equal  amounts  of  the  four  types  of  antiphase  domains  in  the  precipitates.  The 
simulated  evolution  of  alloy  Ni-lOAl-lOCr  at.%  exhibits  the  typical  steady 
nucleation/growth  regime  followed  by  a  coarsening  regime,  where  the  precipi¬ 
tate  area  fraction  approaches  its  equilibrium  value.  The  ordering  profiles  for  Al 
and  Cr  are  shown  separately  in  Fig.  5  at  two  different  stages  of  the  evolution, 
where  the  four  types  of  ordered  domain  structures  having  different  combina¬ 
tions  of  the  three  Iro  parameters  (as  discussed  in  Sec.  2.1.2)  are  represented 
with  four  different  colors  (red,  blue,  white,  aquamarine)  and  the  disordered 
matrix  phase  is  in  black.  The  profiles  of  the  ordered  domains  for  the  two  solute 
elements  overlap  each  other  during  the  entire  evolution  as  shown  in  the  repre¬ 
sentative  snapshots  in  Fig.  5.  These  proof-of-principle  tests  demonstrate  that 
the  model  can  accurately  describe  the  compositional  and  structural  changes 
associated  with  the  precipitation  and  evolution  of  the  7'  phase. 

It  has  long  been  recognized  that  the  ternary  diffusion  behavior  is  quite  differ¬ 
ent  from  that  of  the  binary  systems.  A  classic  example  to  demonstrate  this 
difference  is  the  so-called  uphill  diffusion,  which  represents  a  situation  where 
the  diffusing  substance  can  diffuse  toward  a  higher  concentration  region.  The 
uphill  diffusion  is  widely  observed  in  system  of  more  than  two  components 
since  the  classic  Darken  experiment  in  steels  [24],  In  binary  systems,  the  dif¬ 
fusion  is  driven  by  the  concentration  gradient.  The  species  diffuse  toward  a 
region  of  lower  concentration  and  therefore  uphill  diffusion  is  not  expected. 

In  ternary  or  higher  order  systems,  the  driving  force  for  diffusion  is  attributed 
to  the  negative  gradient  of  the  chemical  potential  of  the  diffusing  substance  [24], 
While  the  concentration  gradient  of  the  diffusing  substance  certainly  con¬ 
tributes  to  the  gradient  of  the  chemical  potential,  the  presence  of  other  al¬ 
loying  components  could  alter  its  chemical  potential  as  well.  As  a  result,  the 
diffusion  path  of  a  species  in  a  ternary  or  higher  order  systems  is  dependent 
on  the  concentration  gradients  of  all  the  constituent  components. 

To  examine  the  capability  of  the  phase-field  method  in  capturing  the  essential 
physics  of  ternary  diffusion,  we  choose  a  few  diffusion  couples  to  simulate 
the  evolution  of  concentration  profiles  under  isothermal  heat  treatment.  The 
following  three  types  of  diffusion  couples  are  studied  [at.%]: 

(1)  Ni-15Al-lCr/Ni-15Al-5Cr; 

(2)  Ni-15.5Al-lCr/Ni-15Al-5Cr; 

(3)  Ni-18Al-lCr/Ni-15Al-5Cr. 

In  all  these  diffusion  couples,  the  right  end  of  the  alloys  are  the  same,  i.e.  Ni- 
15Al-5Cr.  The  Cr  content  are  held  constant  (1  at.%)  for  the  alloys  on  the  left 
end  while  the  Al  concentration  varies.  In  couple  1,  the  Al  content  is  15  at.% 
that  is  the  same  as  in  the  alloy  on  the  right  end.  Therefore  there  exists  only  a 
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gradient  in  Cr  concentration  but  not  in  A1  across  the  couple  interface.  A  slight 
increase  in  A1  concentration  is  introduced  in  couple  2  leading  to  an  opposite 
sign  in  the  concentration  gradient  of  A1  and  Cr  across  the  couple.  In  couple  3, 
the  concentration  of  A1  is  further  increased  leading  to  a  higher  concentration 
gradient  in  A1  across  the  couple.  At  the  beginning  of  the  simulation,  two 
long  slabs  (each  represented  by  a  200x50  uniform  mesh)  of  alloys  are  ’welded’ 
together  to  form  a  400x50  mesh  for  the  simulation  under  isothermal  holding. 
The  alloys  on  both  ends  of  the  couples  are  homogeneous  in  concentration  and 
have  a  single  phase  of  7  for  the  sake  of  simplicity. 

Fig.  6  shows  a  snapshot  of  the  concentration  of  A1  (the  faster  diffusing  species 
relative  to  Cr  in  Ni-Al-Cr  systems)  as  a  function  of  the  distance  from  the  orig¬ 
inal  interface  after  isothermal  annealing  at  873K.  The  positive  (negative)  sign 
of  the  distance  indicates  the  right  (left)  hand  side  of  couple.  The  concentration 
shown  is  an  average  value  of  all  grid  points  at  a  given  distance  from  the  inter¬ 
face.  Only  the  diffusion  zone  (central  portion  of  the  simulation  cell)  is  shown. 
It  is  interesting  to  observe  uphill  diffusion  in  couple  1  despite  the  absence  of  a 
A1  concentration  gradient.  Similar  uphill  diffusion  has  been  observed  recently 
for  Ru  diffusion  in  a  Ni-Re-Ru  diffusion  couple  (Ni-lORe-lORu/Ni-lORu)  ex¬ 
periment  [25].  They  related  this  uphill  diffusion  of  Ru  to  the  negative  cross 
term  of  the  diffusion  coefficien  (DR^Re),  which  promoted  the  more  mobile  Ru 
atoms  to  diffuse  to  the  Re  rich  side  of  the  interface  [25].  We  will  examine  this 
factor  later.  In  couple  2,  the  A1  concentration  on  the  left  side  alloy  is  slightly 
higher.  If  this  is  a  binary  system  without  Cr  presence  on  both  ends  of  the 
couple,  one  would  expect  the  A1  concentration  goes  lower  across  the  couple 
from  the  left  to  the  right.  In  the  current  situation  with  the  presence  of  Cr  on 
both  ends,  we  observe  the  uphill  diffusion  instead.  With  the  A1  concentration 
further  increases  for  the  alloy  on  the  left  end  in  couple  3,  the  gradient  of  A1 
is  further  increased  across  the  couple.  This  increase  in  the  A1  concentration 
gradient  is  enough  to  overcome  the  effect  of  Cr  concentration  gradient  on  the 
chemical  potential  imbalance  across  the  couple  leading  to  the  normal  diffusion 
of  A1  across  the  couple  interface. 

Fig.  7  presents  the  chemical  potential  profiles  across  the  interface  for  A1  and 
Cr  for  the  three  couples.  As  can  be  seen  for  couple  1  in  Fig.  7(a),  the  chemical 
potential  of  A1  exhibits  a  positive  correlation  with  the  Cr  content.  With  lower 
Cr  concentration  on  the  left,  chemical  potential  is  lower  on  the  left  half  and 
higher  on  the  right  leading  to  the  uphill  diffusion.  While  the  A1  chemical 
potential  gradient  decreases  in  couple  2  due  to  higher  concentration  of  A1  on 
the  left  side  of  the  couple,  the  alloy  on  the  left  still  has  a  higher  chemical 
potential  for  Al.  As  the  A1  concentration  is  further  increased  to  18. Oat. %  for 
the  alloy  on  the  left  side  in  couple  3,  which  represents  a  3%  difference  in  Al 
concentration  across  the  interface  of  the  couple,  the  sign  of  the  Al  chemical 
potential  gradient  is  changed  as  compared  to  those  in  couples  1  and  2.  This 
shows  that  a  3%  difference  in  Al  concentration  has  overcome  the  effect  of  a 
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4%  Cr  concentration  difference  across  the  couple.  As  such,  the  uphill  diffusion 
is  not  observed  in  diffusion  couple  3.  As  for  the  chemical  potential  of  Cr 
(Fig.  7(b)),  the  increase  in  A1  concentration  on  the  left  has  decreased  its 
gradient  across  the  interface  (compare  the  three  couples),  but  the  4%  difference 
in  Cr  concentration  on  both  ends  is  significant  and  therefore  the  chemical 
potential  is  always  lower  on  the  left  side  with  a  lower  Cr  concentration. 


We  now  examine  the  effect  of  the  diffusion  coefficients.  According  to  the  Fick- 
Onsager  law  [26,  27],  the  1-dimensional  (1-D)  diffusion  of  A1  in  the  ternay 
diffusion  couple  can  be  described  as: 
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(28) 


where  Daiai  is  the  interdiffusion  coefficient  describing  the  influence  of  A1  con¬ 
centration  gradient  on  its  own  fluxes,  DAicr  the  cross  term  of  the  inter  diffusion 
coefficients  reflecting  the  influence  of  the  concentration  gradient  of  Cr  on  the 
flux  of  Al.  In  our  phase-field  formulation,  the  diffusion  mobilities  are  used  to 
play  the  role  of  the  diffusion  coefficients  as  shown  below  the  1-D  version  of 
the  diffusion  equation  for  Al  from  Eq.  11: 


dxAi 

dt 


d_ 

dx 


d_ 

dx 


| '  MaiCt 


d(6F/5xCr)\ 

Yx  )\  ’ 


(29) 


where  5F/5xai  and  5F/5xcr  are  the  chemical  potential  for  Al  and  Cr,  respec¬ 
tively.  Similar  to  the  interdiffusion  coefficients,  Maiai  is  the  diffusion  mobility 
that  characterizes  the  influence  of  Al  chemical  potential  gradient  on  its  own 
fluxes,  Maict  the  cross  term  of  the  diffusion  mobility  that  describes  the  influ¬ 
ence  of  Cr  chemical  potential  gradient  on  the  flux  of  Al. 


Fig.  8  shows  the  profiles  of  three  diffusion  mobility  terms,  i.e.  Maiai ,  Maict , 
and  Merer ,  across  the  interface  for  couple  1.  As  can  be  seen  from  the  plot,  the 
value  of  Maiai  is  about  an  order  of  magnitude  larger  than  that  of  Mcrcr  con¬ 
firming  that  Al  is  the  faster  diffuser.  The  cross  term  of  the  diffusion  mobility 
Maict  is  negative  and  its  value  is  more  than  an  order  of  magnitude  smaller 
than  that  of  Maiai ■  The  negative  sign  means  that  the  Cr  chemical  poten¬ 
tial  gradient  (mainly  controlled  by  the  cr  concentration  gradient)  in  couple  1 
would  promote  Al  to  diffuse  from  the  left  to  the  right  side  due  to  a  negative  Cr 
concentration  gradient  across  the  interface.  Therefore,  this  cross  diffusion  is 
not  responsible  for  the  uphill  diffusion  observed  in  couple  1.  Also  the  relatively 
small  value  of  Maict  makes  the  direct  contribution  of  the  Cr  concentration 
gradient  to  Al  diffusion  far  less  significant  than  that  of  its  effect  on  the  chem¬ 
ical  potential  gradient  across  the  interface. 


Here  we  discuss  briefly  some  potential  applications  of  this  model.  We  have 
demonstrated  that  the  model  can  accurately  describe  the  precipitation  of  7' 
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in  the  7  matrix  of  Ni-based  superalloys.  This  model  can  be  applied  to  study 
the  precipitation  behavior  in  Ni-Al-Cr  system  quantitatively.  For  example, 
by  examing  the  evolution  of  the  concentrations  in  the  precipitate  under  well 
designed  simulation  conditions  such  as  various  initial  concentration/ordering 
profiles  for  the  nucleus  of  precipitate  may  provide  some  insight  to  the  so- 
called  Cr-trapping  behavior  observed  experimentally  [28],  a  situation  where 
the  concentration  of  Cr  in  the  precipitate  decreases  monotonically  with  the 
growth  of  the  precipitate  while  the  CALPHAD  approach  predicts  the  opposite 
trend.  Details  will  be  published  in  a  forthcoming  paper.  This  model  can  also 
serve  as  the  basis  for  a  pseudo-ternary  model  that  can  simualte  precipitation 
behavior  in  Ni-based  superalloys.  Also  we  treated  explicitly  the  misfit  effect 
due  to  partitioning  of  A1  and  Cr  in  our  elastic  energy  formualtion.  In  this 
work,  we  set  the  misfit  to  zero  since  no  significant  misfit  is  reported  for  Ni-Al- 
Cr  systems.  But  this  formulation  is  general  and  could  be  useful  for  studying 
the  effect  of  misfit  arising  from  multiple  sources  for  alloy  designing  purpose. 
We  studied  the  uphill  diffusion  without  involving  7'  microstructure.  But  it  is 
straightforward  to  incorporate  a  dynamically  evolving  7  and  7'  microstrucure 
on  both  side  of  the  diffusion  couple  and  study  how  different  microstructure 
characteristics  (such  as  phase  fraction  and  their  size  distribution)  affect  the 
diffusion  path  [29].  Considering  that  we  have  implemented  the  most  up-to- 
date  mobility  database,  the  accuracy  of  diffusion  path  prediction  might  be 
improved. 


4  Conclusions 


A  ternary  phase-field  model  based  on  Ni-Al-Cr  systems  was  developed.  The 
model  is  linked  directly  to  commercial  CALPHAD  software  and  implemented 
a  diffusion  mobiltiy  database  that  became  available  recently  including  the  con¬ 
tribution  from  the  ordered  phase.  The  direct  linkage  with  commercial  CAL¬ 
PHAD  software  provides  quantitative  thermodynamic  driving  forces.  And  a 
more  accurate  description  of  the  diffusion  behavior  is  achieved  in  the  alloy  sys¬ 
tem  by  including  contribution  from  the  ordered  state.  Simulations  show  that 
the  model  can  accurately  represent  the  compositional  and  structural  changes 
associated  with  the  7  — >-  7'  transformation  in  the  Ni-based  superalloys.  Uphill 
diffusion  of  A1  is  observed  in  a  few  diffusion  couples  showing  the  significant 
influence  of  Cr  on  the  chemical  potential  of  Al. 

This  model  can  be  applied  to  study  the  precipitation  behavior  in  Ni-Al-Cr  sys¬ 
tem  quantitatively.  This  model  can  also  serve  as  the  basis  for  a  pseudo-ternary 
model  that  can  simulate  precipitation  behavior  in  Ni-based  superalloys  includ¬ 
ing  the  misfit  effect  due  to  partititioning  of  Al  and  Cr  using  the  elastic  energy 
formualtion.  Given  the  most  up-to-date  mobility  database  that  is  incorporated 
in  the  model,  this  model  might  be  employed  to  study  the  microstructure  effect 
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on  diffusion  path  with  higher  accuracy. 
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Fig.  3.  Growth  of  the  precipitate  from  a  small  seed  and  the  accompaning  evolution  of 
concentrations  and  a  representative  Iro  parameter  r/1  for  Al  and  Cr  in  Ni-lOAl-lOCr 
at.%  alloy.  Note  different  y-scale  is  applied  in  (c)  for  the  ordering  and  concentration 
profiles. 


Fig.  4.  Evolution  of  concentration  and  Iro  parameter  r/1  for  Al  and  Cr  in 
Ni-15Al-5Cr  at.%  alloy  at  two  different  stages  of  the  growth  of  a  precipiate. 
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Fig.  5.  Snapshots  of  the  ordering  domains  at  t= 10,  50,  150  (hrs.)  (left  to  right) 
during  the  7'  precipitation  from  the  7  matrix  in  Ni-lOAl-lOCr  at.%  alloy.  Top  row 
shows  the  ordering  profiles  of  Al  and  the  bottom  Cr. 
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Fig.  7.  Chemical  potentials  of  A1  (a)  and  Cr  (b)  for  the  three  couples  with  the 
concentration  profile  shown  in  Fig.  6.  Note  that  the  scale  in  the  two  figures  are 
different. 
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Fig.  8.  Profiles  of  diffusion  mobilities  for  couple  1  corresponding  to  the  concentration 
profile  shown  in  Fig.  6 
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